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A Theta lift representation for the 
Kawazumi-Zhang and Faltings invariants of 
genus-two Riemann surfaces 


By Boris Pioline , Geneva 


Abstract. The Kawazumi-Zhang invariant p for compact genus-two Riemann surfaces 
was recently shown to be an eigenmode of the Laplacian on the Siegel upper half-plane, away 
from the separating degeneration divisor. Using this fact and the known behavior of ip in the 
non-separating degeneration limit, it is shown that p is equal to the Theta lift of the unique (up 
to normalization) weak Jacobi form of weight —2. This identification provides the complete 
Fourier-Jacobi expansion of p near the non-separating node, gives full control on the asymp¬ 
totics of p in the various degeneration limits, and provides an efficient numerical procedure 
to evaluate p to arbitrary accuracy. It also reveals a mock-type holomorphic Siegel modular 
form of weight —2 underlying p. From the general relation between the Faltings invariant, the 
Kawazumi-Zhang invariant and the discriminant for hyperelliptic Riemann surfaces, a Theta 
lift representation for the Faltings invariant in genus two readily follows. 


1. Introduction 


The Kawazumi-Zhang invariant, introduced in [25,32], is areal-valued function </?(£) on 


the moduli space M. h of compact Riemann surfaces £ of genus h > I. One way of defining it 


is through the spectrum of the Laplacian As with respect to the Arakelov metric on £, 



( 1 ) 


where (wi,..., Uh) is an orthonormal basis of holomorphic differentials on £, 0 = Ao < Ai < 
A 2 < ... are the eigenvalues of As, and (f>£ a corresponding orthonormal basis of real square- 
integrable eigenmodes. For genus 2, y>(£) = p(ft) is a function of the period matrix 0, 
and defines a real-analytic modular function on the Siegel upper half-plane 'H 2 , away from 
the separating degeneration divisor. The Kawazumi-Zhang invariant is a close cousin [11] 
of the Faltings invariant <5^(£) [17], which plays an important role in arithmetic geometry. 
Its asymptotic behavior near the boundaries of the moduli space Mh in arbitrary genus was 
investigated in [10,12,30]. 

While Faltings’ invariant made an appearance in studies of bosonisation in conformal 
field theory [1], the genus-two Kawazumi-Zhang invariant has entered the physics literature 
in a recent analysis of the low energy expansion of the two-loop four-graviton amplitude in 
superstring theory [13]: the leading ZJ 4 72 4 interaction is proportional to the Weil-Petersson 
volume of the moduli space M 2 , while the next-to-leading D 6 72 4 interaction is proportional to 
the integral of p times the same Weil-Petersson volume form on M. 2 - With hindsight from var¬ 
ious physics conjectures, it was proven in [14] that p is an eigenmode of the Laplacian Asp(4) 
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on the Siegel upper half-plane, up to a source term supported on the separating degeneration 
divisor, 

(2) [A 5 p( 4 ) — 5] <p = —2ir det(Imfi) , 

where v is the off-diagonal element in the period matrix 


(3) 






As we shall see, this partial differential equation provides strong constraints on the asymptotic 
behavior at the boundaries of A4 2 - 

In the maximal degeneration limit, where all entries in the imaginary part < >2 of the period 
matrix fi = fii + 1 Q 2 are scaled to infinity, it was shown in [10,14] that 


(4) 




7 r 
6 


L\ + L 2 + 


_5L1L2-L3_ 

L1L2 + L2L3 + L3L1 


+ 0(1 /Li) 


where 0 < L 3 < L\ < L 2 parametrize the imaginary part of Q in the standard fundamental 
domain of the action of GL( 2, Z) on the space of 2 x 2 positive definite real matrices, 


(5) 


{Li + L3 L3 \ 

\ L 3 L 2 + J 


This parametrization is motivated by the connection to two-loop supergravity amplitudes, where 
the Li s play the role of Schwinger time parameters [14,19,20], The leading term in (4) is an 
exact solution of ( 2 ) with no source term, which was one of the hints towards the differential 
equation (2) in [14], 

In the minimal degeneration limit a —>• ioo, keeping the other entries of 12 fixed, one has 
instead [12,30] 

( 6 ) p(Q) = + (p 0 (p,u 1 ,u 2 ) + (D(l/t) 

0 

where 

0(p, PU 2 ~ to) ' 

v(p) 

Here, v = pu2 — u\ where ui,U2 are real, t = 02 — is non-negative, r\(p) is Dedekind's 

eta function and 9 (p,v ) = Xmez p+ 2 m{n+ 1 )(v+- 2 ) - s j aco ^p s theta series. The first 

two terms in ( 6 ) satisfy (2) up to terms of order 1/t. Up to the order displayed, each term in 
the Laurent expansion around t = 00 is a real-analytic function of p. u\. 11,2 invariant under the 
Jacobi subgroup Tj = SL( 2, Z) x Z 2 x Z of the Siegel modular group T = Sp( 4, Z) (i. e. a 
Jacobi form of weight zero and index zero). As we shall see, this structure extends to all orders 
in 1 ft. 

Finally, in the separating degeneration limit v —> 0, keeping p, a fixed, one has [12,30] 

<p(ft) = ~ log |27 tz; ? 7 2 (/o)77 2 (cr)| + 0(\v\ 2 log |n|) . 


(V) 


(p 0 (p,U 1 ,U 2 ) = - log 


e -7rp 2 «2 


(8) 
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Eq. (8) is consistent with the differential equation (2) at the order stated, with the logarithmic 
behavior at v = 0 being responsible for the delta-function source term. In [14] it was shown 
using the properties above that the average value of p on M 2 with respect to the Weil-Petersson 
volume form is equal to 3/2, verifying a prediction from S-duality in superstring theory. 

Our goal in this work is to determine the complete asymptotics of the invariant p(Q) in 
the degeneration limits (4) and (6), and more generally, obtain the complete Fourier expansion 
with respect to Oi. To this aim, in §3.2 we shall construct a real-analytic Siegel modular form 
</$(0) on M .2 that satisfies (2), (4) and (6). Since p — p is square integrable on M 2 and an 
eigenmode of ^Sp( 4 ) with non-negative eigenvalue, it must therefore vanish (Theorem 1 in 
§3.3). p is constructed as the Theta lift of the unique weight —2 weak Jacobi form d 2 /tf (see 
Eq. (41)). This parallels the construction of the log-norm of the Igusa cusp form T'io as the 
Theta lift of the unique weight 0 weak Jacobi form (also known as the elliptic genus of K3) due 
to Kawai [23], which we review in §3.1. Since singular Theta lifts were studied extensively 
in [7,9], we refer to these works for issues of convergence. For the convenience of the reader 
however, we shall rederive the Fourier expansions at a physicist’s level of rigor. Using the 
relation between the Fairings invariant, the Kawazumi-Zhang invariant and T 10 established in 
[11], a Theta lift representation for the Fairings invariant is readily obtained (Corollary 3 in 
§3.3). 

The Theta lift representation of p has several interesting consequences, considered in 
§3.3 and §4. First, it gives complete control over the asymptotics in the various degeneration 
limits, and provides an efficient numerical procedure to evaluate p to arbitrary accuracy. This 
is likely to have useful applications in Arakelov geometry. Second, it reveals a ‘holomorphic 
prepotential’ F\ (fl) which generates p through the action of (the real part of) the Siegel-Maass 
raising operator (Eq. (83)). F\ transforms non-homogeneously under Sp( 4, Z), giving an 
explicit example of a mock-type Siegel modular form. Third, it implies that p is an eigenmode 
of an invariant quartic differential operator (Eq. (91)). It would be interesting to understand 
whether the differential equations (2) and (91) can be generalized to higher genus. 

From the physics point of view, the results obtained here will be key for checking S- 
duality predictions for D G 7l l couplings in string theory [29]. In a different vein, it is worth 
noting that the same type of prepotential F\ appears in the physics literature when comput¬ 
ing one-loop corrections to the holomorphic prepotential in heterotic vacua with N = 2 
supersymmetry [5,22,28]. In that context, F\ encodes a subset of the Gromov-Witten in¬ 
variants in the dual type IIA string theory compactified on a suitable K3-fibered Calabi-Yau 
threefold. This analogy suggests that the product of the moduli space of genus-two Riemann 
surfaces times the Poincare upper half-plane FLi (parametrizing the size s of the base of the 
K3-fibration) may cany some canonical special Kahler metric derived from a prepotential 
F(s, p,v,a) = s(pa — v 2 ') + F\ + 0(e~ s ), where F\ is the holomorphic prepotential un¬ 
derlying the Kawazumi-Zhang invariant. It would be very interesting to find a string theory 
compactification whose moduli space carries this putative metric, and compute the 0(e~ s ) 
corrections using mirror symmetry techniques. 


2. Refined degeneration formulae 

In this section, we shall attempt to improve the accuracy of the asymptotic expansions 
(4), (6) and (8) by requiring consistency with the Faplace equation (2) and invariance under 
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the Jacobi group Tj. This section is heuristic, and the proof that p actually satisfies these 
improved asymptotic expansions is deferred to §3. This attempt is inspired by a study of two- 
loop amplitudes in superstring theory [29], and in turn, in combination with insights gained 
from a study of generalized Borcherds lifts [2-4], inspired the educated guess considered in 
§3.2. The reader uninterested by the source of this guess can safely skip to §3. 

Starting with the minimal non-separating degeneration, we observe that the expansion 
( 6 ) can be strengthened, consistently with the Laplace equation (2) to exponential accuracy, to 

(9) = + ~ + 0{e *) , 

b t 

where p± is a function of p, u \, u 2 to be determined. Indeed, decomposing the Laplace operator 
into 


(10) ^Sp( 4 ) — Aj + A p + t A u 4 - A cri 

where 

A t =t 2 d$-tdt, A p = P l(d 2 + d 2 ), A u = ^—\pd Ul +d U2 \ 2 

( 11 ) Z P 2 

A ffl = - (t + p 2 ul)d 2 ai + ( 2tp 2 u 2 d Ul - 2 pluld Pl )d ai , 

and using the fact that pa, defined in (7), satisfies 

(12) A p p 0 = 0 , A„ tp 0 = 7T, 

we see that the Laplace equation (2) is satisfied at order 0{e ~ t ) provided p\ satisfies 


(13) 


Ap pi = 2ipi , A u pi = 5po . 


Invariance of p under T requires that p\ be a real-analytic Jacobi form of zero weight and zero 
index. On the other hand, the maximal degeneration limit (4) requires that, in the limit p 2 —>• 00 
keeping pi,ui,U 2 fixed (with U 2 € [ 0 , 1 ]), 

(14) 990 ~ 7 ^ 2(1 - 6 m 2 + 6 ^) , pi ~ ^rplul(u 2 - l ) 2 • 

fa fa 

The first equation is of course satisfied by (7). It is suggestive to rewrite these limits in terms 
of the Bernoulli polynomials B 2 (x) = x 2 — x + B±(x) = x 2 (x — l ) 2 — 

(15) pa - TTp 2 B 2 {u2) , pi ~ -^-pl {B 4 (u 2 ) + ^j) 

A solution to (13) obeying these boundary conditions can be obtained as a linear combination 


(16) 


v ' = ^p 1 n ^ v) + h E ' {2 ' ll) 


of the standard non-holomorphic Eisenstein series 


(17) 


E‘{s-,p) = \i,-’T{ S ) Y. 

(m,n)p( 0 , 0 ) 


P 2 


| mp + n | 2 
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and the Kronecker-Eisenstein series introduced in [31] 

(2l0^ a+b ~ 1 e 2TT\(nu2+mu 1 ) 

(18) V afi (p-v)= — E (.mp + n) a (mp + n) b ’ 

where a, 6 are non-negative integers. I+j/p: v) is a real-analytic Jacobi modular form of 
weight (1 — b, 1 — a) and zero index, with Fourier expansion 

(19) 

°° °° (4irn 1 

VaM v) = J2 x ) + C- 1 )^ 6 E D ^(Q m X- 1 ) + (! h s, B a+b ( u 2 ) , 

( a + 6 ) ! 

where x = e 2- '™ = e 27n (“ 2 P-ui/ g _ e 2l7r P ? B a (x) are the Bernoulli polynomials, and D ab {x) 

are the Bloch-Wigner-Ramakrishnan single-valued polylogarithms [31], 


( 20 ) 


a+6— 1 

A++) =(-l)“" 1 ^ 2 a+fe - 1 - fc 

k=a 
a+6—1 

I ^_-q6—1 ^ ^ c^a+b—l—k 

k=b 


k — l\ (— log |x| 


\a-\-b—l—k 


■ Lifc(x) 


a—1/ (a + 6—1 — A:)! 


k — 1 \ (— log |x|) a+b 1 k — 


6—1/ (a + 6 — 1 — k)\ 


Li k{x) ■ 


It is easy to check that the differential equations (13) are obeyed, by checking the action on the 
seed of the Poincare series (i.e. setting m = 0, n = 1) and using the second Kronecker limit 
formula, which states 

(21) +o = ^T> h i(p-v) . 

Moreover, the equality (16) predicts an additional subleading term in (15), 


( 22 ) = y/^ 4 (+) + |j Pi + P2 l + 0(e~ 2 ^) 

The third term in (22) requires a subleading correction to the maximal degeneration limit (4), 

(23) 


<p(ty = 


7T 


+1 + +2 + +3 — 


5 + 1 + 2+3 


+ 1+2 + +2L3 + +3 + 1 


+ 


5C(3) 


47T 2 (+i+2 + +2+3 + +3+1) 


+ ... 


The additional term is an exact solution of the Laplace equation (2). 

With the hindsight gained from a study of generalized Borcherds lifts [2], the estimates 
(9) and (23), if true, strongly suggest that tp(Q) is the Theta lift of an almost, weakly holomor- 
phic Jacobi form of weight —1/2 and depth 1, motivating the educated guess in §3.2. We shall 
prove in §3.3 that this estimates do in fact hold with exponential accuracy. 


3. The Kawazumi-Zhang invariant as a Theta lift 

In this section, using a suitable Theta lift, we construct a real-analytic Siegel modular 
form (p that satisfies the differential equation (2) and asymptotic behaviors (4), ( 6 ), ( 8 ) in the 
various degeneration limits - hence must coincide with p. As a warm-up, we start by recalling 
the Theta lift representation of the log-norm of the discriminant of genus two Riemann surfaces, 
following [23]. 
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3.1. The Igusa cusp form T'io as a Theta lift. Recall that the log-norm log |'I ; io| 
= log[(det 0 2 ) 5 |^io|] of the weight 10 Igusa cusp form (normalized as 2 ~ 12 times the product 
of the squares of the ten Thetanullwerte) can be represented as a regularized modular integral ^ 


logH^ioll(O) 

(24) 


1 


L 


d 2 r 


4 

5 log 


Tl r 2 
87T 




3\/3 


a t-7 E 


20 r 2 


where T\ = {r € Tii, |r| > 1 , — -7 < t\ < 4} is the standard fundamental domain for the 
action of SL(2, Z) on "Hi, d 2 r = dridr 2 , 7 e is the Euler-Mascheroni constant, and p® v ® n l odd 
are partition functions for even (shifted) lattices of signature (3,2), or Siegel-Narain theta series, 


y-ievenlodd , \ 

r 3,2 (^;r) =t 2 2 ^ 


in 2 -in 2 


(25) 


(mi,m 2 ,n 1 ,n 2 )£Z 4 

fee 2Z|2Z+1 

2 |m 2 — p mi + an 1 + (pa — v 2 ) n 2 — 6 n | 2 

PR 2 : 

P2 <?2 ~ 

pi =Pr + 4 min 1 + 6 2 . 


For v —» 0, T| v | n —> r 2 , 2 0 3 (2r),r§ d 2 d -A T 2 ,2 02(2t) where T 2 ,2 is the partition function of 
the even self-dual lattice with signature ( 2 , 2 ), 


(26) 


r 2 ,2 (p,ct;t) = t 2 ^ 

(mi ,m2,n 1 ,n 2 )GZ 4 


I „2 _i 2 

g 4 Pi g4Pfl | 6=1)=0 


Thus, p®™ n l odd are modular forms of weight 1/2 under To(4). Under general modular trans¬ 
formations of r, 


r|3”(r +1) =rjj”(T), rg>(r + 1 ) = i r^(r) 


~iodd / 


(27) 


1 -i 


r|;r(~i/r) =— T 1 / 2 [r-| n (r) + r° dd (r)] , 


r 3 d 2 d (-l/ r ) = V rV2 FITM - r 3 d 2 d (^)] • 


The integers mi, m 2 , n 1 , n 2 , 6 can be fit into an antisymmetric traceless matrix 


(28) 


0 

—m 2 

6/2 

n 1 ^ 

m 2 

0 

mi 

-6/2 

—6/2 

—mi 

0 

—n 2 

—n 1 

6/2 

n 2 

0 / 


with Pfaffian proportional to p 2 L — pf { , which transforms by conjugation under Sp( 4, Z).This 
makes it clear that Sp( 4, Z) transformations preserve the parity of b. Thus, both p;/'/"l' ,dd (Q: r) 

11 The formula (24) was discovered in [23] by computing threshold corrections to gauge couplings in het¬ 
erotic string theory compactified on A'3 x T 2 . The variables p, a, v parametrize the complex structure, Kiihler class 
and holonomies of a (7(1) connection on the torus T 2 , while the integers mi,n t correspond to the momentum and 
winding numbers, and b is the electric charge. 
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are Siegel modular functions in the variable U, and so then is the result of the modular integral 
(24). On the other hand, ho, hi are the coefficients of the theta series decomposition of the 
elliptic genus of K 3, 


(29) 

XK 3 (r, z) =h Q {r) 9 3 (2 t, 2 z) + hi(r) 9 2 (2 t, 2 z) , 

ho(r) =24 - 2 ^ r) ~ 9 2 (2t) = 20 + 216 q + 1616 q 2 + ... 


9 K r ) 


V e (r) 


hi(r) = 24 ^^ + 2 ^ (t) 6 f [t) 9 3 (2t) = q" 1 ' 4 (2 - 128 q - 1026 q 2 + ...) 
0|(r) r/ b (r) 


They are modular forms of To (4) with weight —1/2. In terms of the standard generators 
X 2 {2t) = E 2 (2t) — 2E 2 (At), 9 2 (2t) of the ring of To(4) modular forms of even weight, 


(30) 


ho 


hi 


1 

9 3 {2t)A 6 

1 

9 2 (2t)A 6 


■^9?(2t) - X 2 (2t)91{2t) + ^[X 2 (2t)] 2 0 4 (2t) 

- A 0 12( 2t ) + x 2 (2t)9 8 2 (2t) + \[X 2 {2t)} 2 0 4 (2t) 


where A 6 = r/(2r)] 12 is a cusp form of weight 6. Under general modular transformations, 


h 0 {T + 1) =ho(r) , hi(r + 1) = -i/ii(r) , 

(31) hoi-l/r) =^~y- r~ l/2 Mu) + M T )] , 

M-1/ t ) =—E~ r_1/2 [ h o ( T ) “ M T )] > 

so that the integrand of (24) is (except for the last term in the bracket, proportional to r 2 ) 
invariant under the full modular group SL( 2,Z). It follows from (31) (see e.g. the proof of 
Thm 5.4 in [16]) that 


(32) 

h(r) =/i 0 (4r) + hi (At) 

= - 803(2 1 r)A6 ( — 30 2 2 (2t) + 40|(2t)A 2 (2t) + 120 4 (2t)X 2 (2t) - 16A|(2r)) 

= c(m)q m = 2q- 1 + 20 - 128<? 3 + 216g 4 - 1026q 7 + ... 

m> — 1 

is modular form of weight —4 under To (4) in Kohnen’s plus space (;. e. the m-th Fourier 
coefficient c(m ) of h(r) vanishes unless m = 0, 3 mod 4). h(r) has a simple pole at r = ioo 
and is regular at r = 0 and r = /. The constant term in h(r) makes it necessary to subtract by 
hand the term proportional to r 2 in (24), in order for the integral to converge. 2} 

Using the differential equation satisfied by the lattice partition function, 

(33) [As p (4) - 4A5 i ( 2 ) il / 2 + l] r^ n|odd = 0 , 

2) Alternatively, following [7] one could truncate the integration domain to T\ = Ei fl {t 2 < A}, insert 

a Kronecker regulating factor r| in the integrand, take the limit A —r oo for fixed s with Re(s) sufficiently large, 

analytically continue in s and extract the constant term in the Laurent expansion at s = 0. The two prescriptions 

can be shown to agree up to an additive constant. 
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where A sl( 2 ),w = 4r|9^(i9 T + 7 ^) + w is the Laplacian acting on modular forms of weight 
w, one sees that log ||4'io|| is areal-analytic quasi-harmonic function on the Siegel upper half¬ 
plane, up to a delta function source term supported on the separating divisor, 

(34) A Sp(4 ) log ||^io11 = -15 + 47r<5 (2) (n) . 

Indeed, as v -A 0, the integrand becomes 

(35) r 2 ,2 M2t) /t 0 (r) + 9 2 (2t) Mr)] - 20t 2 = 24T 2i2 - 20t 2 T2 ~°° 4t 2 , 


which leads to a logarithmic divergence. Keeping v small but non zero, and retaining the 
contributions from m, = n r = 0,b = ±1, we have 


(36) 


log 11 'b 


101 


i; 


dr 2 - 


TT 2 \vr 


T2 


e p ^- v 2 =-T(0,irz) , 


where z = • Using the fact that the incomplete Gamma function F(0. ttz) behaves as 

— log(/rz)+ analytic as z —>• 0, and the result from [15] 

87re 1_7 

o ~(f 2,2 — T 2 ) = — lOg 

Ti ^ 


(37) 

we find 

(38) 


p2V2\ri(p)ri{v)\ 4 


^ b L 3\/3 

log 11^ 10 11 = log|p 2 cr|n 2 r/ 24 (p)r/ 24 (cr)| + 0(\v\ 2 ) , 


where the omitted terms vanish analytically as v -A 0. 

Evaluating the modular integral by the standard unfolding method [7,22], one arrives at 

log 114'io||(U) = - 2n(p2 + cr 2 — v 2 ) + 51ogdetU 2 
( 39 > - Re [ C ( 4M - ft2 ) lo S (l - e 2n ^ ka+ip+bv ^ ] , 

(k,i,b)> 0 

where (k,£, b) > 0 stands for {(£; > 0, i > 0) or (k > 0, l > 0), b € Z}U{fc = £ = 0, b > 0}, 
and Q is assumed to be such that ka 2 + ip 2 + bv 2 > 0 for all ( k , £,b) >0 [23, Eq. (20)]. Eq. 
(39) is consistent with the Gritsenko-Nikulin product formula [21] 

(40) ^io(J2) = e 2?ri (p+cr-v) ^ _ e 2iri(ko+tp+bv)y(Ake.-b 2 ) _ 

(k,i,b)> 0 


3.2. An educated guess. Motivated by the heuristic considerations in §2, and in anal¬ 
ogy with the Theta lift representation of log ||'Fio|| reviewed in §3.1, we consider the modular 
integral 


i r j 2 

(41) <p(n) = -~ -f r|y| n ((l; r) D T ho(r) + r) D T h\{r) 

2 -hi r 2 


where D T = ^(<9 T — ///) is the raising operator, mapping modular forms of weight w to 
modular forms of weight w + 2, and (h°, h 1 ) are weight w = —5/2 modular forms of To(4), 
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associated to the weak Jacobi form <f> = 9 2 (t. z)/ rf (r) of weight —2 and index 1 in the same 
way as before: 


(42) 


4 >{t, z) =ho(r) 9 3 ( 2 t, 2 z) + hi(r) 9 2 { 2 t, 2 z) 

r , v 6» 2 (2r) o 

h 0 (r) = 2 + 12? + 56 q 2 + ... , 


T 


hi (r) = - 


#3 (2r) 




= -g" 1/4 (l + 8q + 39q 2 + ...). 


The corresponding weight —5/2 modular form in Kohnen’s plus space is 

(43) 


[?7(4r)] 6 40 3 (2 t)A 6 

^ c(m)q m = -(--2 + 8g 3 - 12q 4 + 39<? 7 + ... 


m> — 1 


Q 


As in the previous case, h(r) has a simple pole at the cusp at infinity and is regular at the other 
cusps r = 0 and r = The action of the raising operator D r on hi evaluates to 

(44) D r hi = E 2 hi - hi , i = 0,1, 

where E 2 = E 2 — is the almost holomorphic Eisenstein series of weight 2 and h j is defined 
in (30). The zero-th Fourier coefficient term of D r ho is —5/(27 tt2), so the integral (41) is 
convergent, with no need for regularization. Using (33) and the fact that D r hi is an eigenmode 
of Asl( 2 ),-i /2 with eigenvalue 5/2, one easily checks that (p is an eigenmode of A 5p ( 4 ) with 
eigenvalue 5, away from the separating degeneration divisor v = 0. In the limit v -A 0, the 
integrand becomes 


(45) 


2,2 


0 3 ( 2 t) D t Ji 0 (t) + 9 2 ( 2 t) D t Jii(t) 


^ T2-300 

= —1 2,2 ~ — T 2 


leading to a logarithmic divergence. Keeping v small but non zero, retaining the contributions 
from rrii = n l = 0, b = ±1 and using D r h\ ~ — 2g _1 / 4 (l — qAE-), we have 


(46) 



dr 2 - 

-e 

T~2 


71-T2M 2 

--T 

P 2 C7 2~ v 2 



(! + ^)T(0,ttz) - 2 i - e “ 


where z = ——— 


Thus, in the separating degeneration, we have, in agreement with (8), 


(47) <p{Sl) 


l( 1+ 5 H 2 ^ 

2 V ' 2(p 2 (j 2 — v 2 )) 


log |n| 2 


log \2nri 2 (p)rj 2 (a)\ + 0(\v\ 2 ) , 


up to terms vanishing analytically as v —> 0 (the last term on the right-hand side follows from 
(37)). The logarithmic singularity implies that 


(48) [ A S P ( 4 ) - 5] <p = -27t det(U 2 ) 5 (2) (n) , 


so (p satisfies the same equation (2) as <p. In the remainder of this subsection we extract the 
asymptotics of <p in the minimal and maximal non-separating degenerations, and find that they 
agree with the asymptotics of tp. 
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Maximal non-separating degeneration. The maximal degeneration —> oo corre¬ 
sponds, in string theory parlance, to the limit where one of the circles in the torus T 2 becomes 
infinitely large (see footnote 1). In this limit, the lattice partition function r]//"^ Kl ‘ i (Q) factor¬ 
izes into Ti i(r; r)x P'//"' 0 ' 1 '* (f; r), where r = \/det Qo parametrizes the radius of the large 
circle, and f = u 2 + \\Jtj p 2 = Y + \R the radius R and Wilson line Y for the circle of finite 
size. It is useful to express the lattice partition functions Ti i and F^.i in the ‘Lagrangian’ 
representation, where modular invariance in r is manifest, 

r M (r;r)=r 5^ e ~r 2 \p+qr\ 2 /r 2 

(p>9)GZ 2 

(49) -peven|odd/~ , _p —ttR? l m + nT l +2mn(m+nT)Y 2 + i:! ^ L b 2 +2iTr(m+nT)bY 

1 2,1 ( r ; T ) 2-j e 2 

(m,n)sZ 2 

be 2Z|2Z+1 


In the limit r —>• 00 , the 0(3, 2, Z) = Sp( 4, Z) symmetry is broken to 0(2,1, Z) = GL( 2, Z), 
acting on the modulus f € 'H 1 by fractional linear transformations, along with the anti- 
holomorphic involution t 1 —> —f. The leading term in this limit originates from the term 

(p,q) = (0,0) in T 1) i(r), 


(50) 


VL = 




ri;r(r; t) D T h 0 (r ) + T^ d 1 d (f; r) D T hi(r) 


To compute this integral, we decompose the sum over (m, n) in T 2 .i into orbits of SL( 2, Z), 
obtaining pj j = 

(m,n) = (0,0), 


obtaining pj j = ■ The first term corresponds to the contribution of the zero orbit 


(51) 


( 0 ) rR f d 2 r 

^ = "T k t 


6» 3 (2r) D t Ji 0 (t ) + 9 2 (2t) D T hi{r) 


nrR 

~ 6 ~ 


since, as already noted in (45), the term in square bracket reduces to —1. The remaining orbits 
(m, 0) with m / 0 contribute 
(52) 

ri? f d 2 T 


V { L = 


2 Js r| 


E' 

ra /0 


J —7rR 2 m‘ 2 /T2 


63 ( 2 t, 2mY) D T ho(r ) + 0 2 (2t, 2mY) D T h\{r) 


where 5 is the strip [—1/2,1/2] x M + . The integral over n picks up the constant term in 
D T ho(T) (corresponding to b = 0) and the polar term in I) T h] (r) (coiTesponding to b = ±1). 
Thus we have 


(53) 


rR rdra v 

2 Jo r 2 


0 —irR 2 m 2 /t 2 


_5_ 

12 

00 


1 - 


^ m ^0 

3 


7TT2 


(2 — 1 


27riray _ ^—2nimY 


) - 24 t 2 ° + 2 


0 27rimY 


2 e 


—27rimy ^ 


-E 

2 ^ 

m=l 


_cos(2nmY) — 1 2cos(27r mb) 


7r 3 i? 3 m 4 


irRm 2 
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Using the identity Lifc(e 27nx ) + (—l) fc Lifc(e 27nx ) = — Bk(x ) for the polylogarithm, 
valid for k integer, 0 < Re(x) < 1, we arrive, in the region — ^ <Y <\,R 2 + Y 2 > 1, at 


(54) ^<4°> + rf) = y 


R + 5 u 2 (|y|-i) 2 + i_- 6|y | + 6 Y 2 


R 3 


R 


Setting p 2 = r/R,U 2 = Y,t = rR, this reproduces the desired behavior (14) in the maximal 
separating degeneration limit ! The square bracket in (54) is recognized as the local modular 
function A(t ) in the two-loop supergravity computation [20, Eq. (3.8)]. 

The subleading terms in the maximal non-separating degeneration limit are obtained by 
restricting the sum Ti i(r) to the orbit representatives (p, 0) with p / 0, and unfolding on the 
strip: 

(55) 

T1 T(R, Y) D T h 0 (r) + T° 2 f(R, Y) D^t) 

p^o 


p-p L = 


L-iY 


9 9 / 

—7rr^p^/r2 


Up to terms of order e r , one can replace T| v f n —>• y/vi, T^ —>• 0 and D r ho(T) by its constant 
term — 5/(27 tt 2), leading to the next-to-leading correction 



consistently with (22). The exponentially suppressed contributions to <p — (Pl — pnl will be 
obtained in the analysis of the minimal non-separating degeneration, to which we now turn. 


Minimal non-separating degeneration. The limit <72 —>• oo keeping other entries of 
Q fixed corresponds, in string theory parlance, to the limit where the volume of the torus T 2 
becomes infinite, keeping the complex structure p and holonomy v fixed (see footnote 1). The 
Siegel modular group T is now broken to the Jacobi subgroup Tj. Following [7,22], the 
Fourier-Jacobi coefficients (i.e. the Fourier coefficients with respect to or) can be extracted by 
applying the orbit method to the lattice partition function written in the ‘Fagrangian’ represen¬ 
tation, obtained from (25) by Poisson resummation in mi, m 2 , 

(57) r 3 e J n|odd (U;r) = t ^ e ^-nG(A) ^ 

Aei? x2 
b£ 2Z|2Z+1 

where 

(58) 

G(A) = —|- 2icrdet A + —b ■ (v A — v A) + — (v 2 A — v 2 A) — 2i^|(ni + 712 p) A 

P2 t 2 P2 P2 Pi 

and 
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The integer matrix A transforms linearly under SL( 2,Z), and belongs to one of three 
different types of orbits. The orbit A = 0 produces, as in (51), 


(60) 




d 2 i 


ITi 


9 3 (2t) D t Jiq(t) + 6 2 (2t) D t Jii(t) 


nt 

T' 


The degenerate orbits with rk.4 = 1 give instead, restricting to the orbit representatives 
(ni, 77 , 2 ) = ( 0 , 0 ), (mi, m 2 ) / ( 0 , 0 ) and unfolding on the strip. 


( 61 ) 





E ■- 

(mi,m2)/0 


7rt|m^ +m2Pp 

T 2P2 


X 


6 3 (2t, 27771772 + 27772771) D r h() + 0 2 (2t , 2mi7t2 + 2?77277 i) D r h\ 


As in (53), the integral over n picks up the constant term in I) T ho (corresponding to b = 0) 
and the polar term in D T h\ (corresponding to b = ±1), leading to 


(62) 

T 2 


_ E 


,p\(cos[ 2 Tr{m\u 2 + 777 , 2771 )] — 1) 2 p 2 cos[27r(mi7i2 + m27ti)] 


= 2 P hi (P> v )+t 


-1 


7T 3 f 2 |777l + m 2 p\ A 


7rf|777i + m 2 p \ 2 


Combining (60) and (62), we reproduce the desired behavior (9) in the minimal degeneration 
limit, with tpo and ip 1 given in ( 21 ) and (16) ! 

For the non-degenerate orbits with det 4/0, the integral can be unfolded on (a double 
cover of) the upper half-plane, at the expense of restricting the sum to 772 = 0 < mi < ri\ 7 M). 
Substituting the Fourier series of ho and hi, the integral over n is Gaussian, while the integral 
over r 2 is of Bessel type. After some algebra (see e.g. [22, A.2]), we arrive at 


(63) 


c(Ak£ — b 2 ) 

k> 0,f>0 

bez 


5 

16ir 2 tp 2 


E>2,2 (x) + ^(4 kl 


b 2 )D h i (x) , 


where x = e 2 m (ka+ep+bv), ant | j) n h ( x ^ are the single-valued polylogaiithms defined in ( 20 ), 

(64) Di } i(x) = 2Re[Lii(a;)] , D 2:2 (x) = -4Re[Li 3 (x) — log |x| Li 2 (x)] . 

The formula (63) holds in the chamber where ka 2 + lp 2 + bv 2 > 0 for all stated values 
of ( k,£,b ). The sum converges absolutely in a neighborhood of the zero-dimensional cusp 
12 = ioo by the same arguments as in [7,9]. 


3.3. The Kawazumi-Zhang and Faltings invariants as Theta lifts. Using the results 
in §3.2, §3.2 and the facts summarized in the introduction, we see that tp = ip — (p is annihi¬ 
lated by A Sp( 4 ) ~ ^ and vanishes up to order 0 (l/t) in the non-separating degeneration limit 
t —>• 00, and up to order O (j v \ 2 log |t;|) near the separating divisor. On the truncated funda¬ 
mental domain T 2 = T 2 FI {t < A, |t;| > 1/A}, where T 2 is the standard fundamental domain 
from [18], one has 


(65) 


/ ( F 2 (* 1 ) = F / 0Asp(4)0 = -R / dip-k dip[ ip k dip , 

Jrt 5 77-/ 5 Jjr A 5 JgjrA 
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where * denotes the Hodge star on TI 2 , and OF^ the boundary of F^- By the above estimates, 
the boundary term vanishes in the limit A —> 00 , while the first term converges to a finite, non¬ 
positive value. Since the left-hand side is non-negative, it follows that (p must vanish. Thus, 
we have shown the 


Theorem 1. The Kawazumi-Zhang invariant p(0.) for compact genus-two Riemann 
surfaces admits the Theta lift representation 


( 66 ) 


T>{^) = - j 


d 2 r 


'F 


rf| n (T!; r) D r h 0 (r) + r^ d 2 d (H; r) D r h^{r) 


where rg V 2 n l odd (n; r) are the Siegel-Narain theta series defined in (25), and (ho,h\) is the 
weight — | vector-valued modular form appearing in the theta series decomposition (42) of the 
weak Jacobi form 6 |2 (r, z)/r/ Q of weight 2 and index 1. 


Corollary 1. p(f-l) satisfies the following improved asymptotics: in the minimal non- 
separating degeneration t —>• + 00 , 

(67) = -xt + To + ~—h @(e t ) , 

6 t 

where po and p\ are defined in (16); in the maximal non-separating degeneration Li —>• + 00 , 


( 68 ) 


7r 
6 


L\ -f L 2 + £3 


_5 L 1 L 2 -L 3 _ 

L1L2 + L2L3 + L3L1 


5C(3) 

47r 2 det FI 2 


+ 0{e~ Li ) ; 


in the separating degeneration v —> 0, 

(69) = -i ^1 + log |r;| 2 - log |27T77 2 (p)?7 2 (cr)| + 0(\v\ 2 ) . 


Corollary 2. p(fl) admits the Fourier expansion 

rot T , I It 5 tt \v 2 \{p 2 - M )(°2 - M) , 5C(3) 

v(W + *2 - M) - T - 


(70) 


167T 2 det FI 2 


+ i ^ (4 k£ - b 2 ) c(4k£ - b 2 ) £>i,i ( ' e 2 ^+e P +bv )) 


^ c(4M - 6 2 ) £> 2,2 (e 2 ^+F+ bv ) ) 


{k,t,b)> 0 


(k,e,b)>o 


where ( k,£,b ) > 0 was defined below (39), and D 1,1 (t), £>2,2 (x) are given in (64). The 
Fourier expansion is absolutely convergent in a neighborhood of the zero-dimensional cusp 
= ioo. 


In [11], a relation between the Kawazumi-Zhang invariant <£>(£), the Faltings invariant 
<5(£) and the discriminant A(£) for hyperelliptic compact Riemann surfaces £ was obtained. 
At genus two, all compact Riemann surfaces are hyperelliptic, and the discriminant is propor¬ 
tional to the Igusa cusp form 'T 10 - Corollory 1.8 in [11] states 

(71) <p(fi) = —3 log H'I'ioll(fi) - - 401og27r . 

Using (24), (66) and (71), we obtain 
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Corollary 3. The Faltings invariant admits the Theta lift representation 


(72) 


s F (n) = 



2E2I11 + 7 h\ 


The Fourier expansion of 6 F (Q) is easily obtained by combining (39), (70) and (71). 


4. Miscellany 


4.1. Numerical applications. The formulae (39) and (70) provide an efficient numer¬ 
ical procedure for evaluating the Faltings and Kawazumi-Zhang invariants to arbitrary preci¬ 
sion. As an illustration, for the curve y 2 + y = x 5 considered in [8], with automorphism group 
Z 5 x Z 2 and period matrix 



(73) 


we find, truncating the sum at k, £,\b\ < 15, 

(74) ip = 0.53801117620500504861... , S F = -16.6790574451477760445 ... , 

where all displayed digits appear to be stable upon varying the truncation. This is consistent 
with the value 5 F = —16, 679 ... which follows from the numerical computations in [8, §4.5]. 

For another example, consider the curve y 2 = x 6 — 1, with automorphism group Dq x Z 2 
and period matrix 



(75) 


Using the same truncation, we find 

(76) ip = 0.59291015631443383207... , S F = -16.3412295821338262636 ... 

Finally, consider the Burnside curve y 2 = x(x 4 — 1), with automorphism group S 4 x Z 2 and 
period matrix 


(77) 



Using the same truncation, we find 


(78) p = 0.51986038541995901150 ... , d F = -16.8264632650009721134 ... 







14 


Pioline, A Theta lift representation for the Kawazumi-Zhang invariant 


4.2. Relation to Gromov-Witten invariants. We note that the Fourier expansion (70) 
is similar to [22, Eq. (A.44)], where modular integrals of the form fj. ^rF$ t + 2 , 2 E 2 F(r ) were 
considered. Here, Vxt+ 2, 2 is a partition of an even self-dual lattice of signature (8 1 + 2,2) and 
F(t) is a weakly holomorphic modular form of weight —At—2. In this context, the analogue of 
the coefficients c(4kl — b 2 ) were identified as the BPS invariants (also known as Gopakumar- 
Vafa invariants, and closely related to Gromov-Witten invariants) counting rational curves in a 
suitable K3-fibered Calabi-Yau threefold with hip = 8 t + 3. It is therefore natural to ask if the 
coefficients c(4 k£ — b 2 ) in the Fourier expansion of the Kawazumi-Zhang invariant count ratio¬ 
nal curves in a suitable Calabi-Yau threefold with hi 7 = 4. Two examples of threefolds with 
/lip = 4 were studied in [6,24] (see also [27]). For the example X{1., 1, 2,6,10)_372 in [6,27], 
the rational curves are counted by the weight —2 Jacobi form —{lE/iE^p + 5/Tf ; /T 4 . 1)/(6 r/ 24 ). 
For the example X (2, 2, 3,3,10)_i32 in [24], they are instead counted by the weight —2 Jacobi 
form —2EiE( } \ /rf A . The Jacobi form relevant for the Kawazumi-Zhang invariant is propor¬ 
tional to the difference of these two, 0 2 /rf = (F 4 F 0 . \ — EqE^ i )/(144 rft A ). It is unclear to 
the author whether the fact that it is a weak Jacobi form (i.e. has h{r) = 0(1/q) rather than 
0(l/q A ) as in the cases studied in [6,24,27]) disqualifies it from counting rational curves. 


4.3. Holomorphic prepotential. It is known from [22] that modular integrals of the 

form 


(79) 



-peven 
1 3,2 


(n-r)H 0 (r) + T o 3 f(n-,T)Hi(r) 


c( 0 , 0 )r 2 , 


where h{r) = h 0 (Ar) + /ii(4r) = Yjm>- K 0 <^<! c(m, E)q m /r| is an almost weakly holomor¬ 
phic modular form of weight —1/2 and depth 1 under To (4) in Kohnen’s plus space, can be 
expressed as 


(80) X = Re Fo + Re m_ 2 Fi — c(0, 0) log det Cl 2 , 


where F (l and F\ are holomorphic functions of Cl. Here, is the Siegel-Maass raising opera¬ 
tor 


( 81 ) 




7T 


d 0 d a - \d 2 + 


i(l — 2 w) fw 


4(p 2 <j 2 ~ vl) 


(^7 + cr 2 d a + p 2 d p + v 2 d v ^j 


which maps Siegel modular forms of weight w to modular forms of weight w + 2. Fo is the 
logarithm of a holomorphic Siegel modular form of weight — 2c(0,0). F\, known as the holo¬ 
morphic prepotential, is ambiguous modulo elements in the kernel /C of the operator Re(d_ 2 ). 
The latter includes quadratic polynomials in (p. a. v) with imaginary coefficients, as well as cu¬ 
bic polynomials of the form {pcr — v 2 ){ap + ft a) where a, ft are imaginary. Since the integral 
X is a Siegel modular function, F\ must transform under 7 € Sp{ 4, Z) as 


(82) F 1 |_ 27 (H) = Fi(H) + F 7 (H), 

where F 7 (H) is an element in 1C. Thus, F\ is a mock-type holomorphic Siegel modular form 
of weight —2. For the modular integral (24), the modular form H{t ) = —\h{r) is weakly 
holomorphic therefore F\ vanishes, while Fq = log 'I' 10 (up to an additive constant). For the 
modular integral (41), H{t ) = —| D T h{r) is the modular derivative of a weakly holomorphic 
form, therefore Fo vanishes [2,4], Since, by Theorem 1, <g(Q) is equal to X for this choice of 
H{t), we have the 
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Corollary 4. The Kawazumi-Zhang invariant ip(£l) is equal to the real part of the action 
of the Siegel-Maass raising operator 0_ 2 on the ‘holomorphic prepotential’ F\ (Q), 

(83) <,? = Re (□_ 2 Xi) 

where 

■ 3 

(84) Fi(fi) = E c(4/e£- 6 2 ) Li 3 ()) _ ^_ pcJ (p + CJ _ 2u ) + C(3) . 

(k,e,b)>o 


Proof. Using [26, A.32] 


(85) D n _ 2n U 2n+1 (x) =Y 


n\(n + r — w)\ 2 2n 2r (kt — b 2 ) 

'-f r\(n — r)\(n — w)\ (7rdetS22) f 


l2 \n—r 


L 


r =0 


with n = 2, where x = e 27ri ( fcer +^+H and 


( 86 ) 


L(r)(z) - E 


(r + m)! 


— 1 ' m!(r — m)!(47r) r 
m=0 v M ' 


■ [Irn z] 


Li 


r+m+l 


(r) 


( e 2 ^) 


log X 
27ri 


is related to the Bloch-Wigner-Ramakrishnan polylogarithm D a j, in (20) via [2] 


(87) 


Z7 r+ i ir+ i(x) = 2 Re 


\{-^r , 

r\ L(r) 



one easily checks that the action of Re(D_ 2 ) on the first term of (84) produces the last two 
lines in (70). The action of the same on the polynomial terms in (84) produces the first line in 
(70). □ 

Remark. More generally, integrals of the form (79), where H (r) is an almost weakly holomor¬ 
phic modular form of weight —1/2 and depth n, can be expressed as 


(88) X = Re n_ 2 r F r ~ c(0,0) log det Q 2 , 

r =0 

where F r are holomorphic functions of 0 known as generalized prepotentials, which transform 
as mock-type Siegel modular forms of weight —2 r [2,4,26]. When /7(r) is obtained by acting 
r times with the raising operator D r on a holomorphic modular form h(r) of weight — 2r — 
then all F r vanish except F n [2,4]. 


4.4. Quartic differential equation. Observe that the Narain partition function satisfies, 
in addition to (33), 


(89) 


□ 2 D 0 — 


1 


16-7T 4 


Ac 


SL( 2),1/2 l ^SL( 2),1/2 


1 


peven|odd _ ^ 

1 3.2 — U ’ 


where is the Siegel-Mass lowering operator (formally independent of w) 


2\2 


(90) □ = -TT 2 (P2(72-V 2 ) 


g-f) _I A 2 _ l _ 

u p u cr 4 u v A ( 2 \ 

4(p2 ct2 - V%) 


(a 2 dv + p 2 dp + v 2 dy) 


which maps Siegel modular forms of weight w to Siegel modular forms of weight w — 2. By 
integration by parts, we conclude that 















16 


Pioline, A Theta lift representation for the Kawazumi-Zhang invariant 


Corollary 5. satisfies the quartic differential equation (away from the separating 

degeneration) 



(91) 


It would be interesting to understand the fate of the differential equations (2) and (91) at 
higher genus. 
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